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We propose an adaptive denoising scheme for poorly modeled non-Gaussian features 
in the gravitational wave interferometric data. Preliminary tests on real data show 
encouraging results. 
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1. Motivation 

Several large-scale interferometric gravitational wave detectors will come on-line 
soon, such as L1GO in the U.S., the French/Italian Virgo project, GEO600 the 
German/British interferometer and TAMA in Japan. Gravitational wave detectors 
produce an enormous volume of output. Data analysis techniques will have to be 
developed to optimally extract the weak signature of a gravitational wave from 
these data. Many of the techniques developed so far are based on matched filtering 
and assume stationary Gaussian noise. 

However, the real data stream from the detectors is not expected to satisfy the 
stationary and Gaussian assumptions. This disparity between standard Gaussian 
assumptions and real data characteristics poses a major problem to the direct appli- 
cation of matched filtering techninues in particular for burst sources such as black 
hole binary quasinormal ringings □ or inspiral waveforms i. In fact, the data from 
the Caltech 40 meter proto-type interferometer has the expected broadband noise 
spectrum, but superposed on this are several other noise features: such as long- 
term sinusoidal disturbances coming from suspensions and electric main harmonics 
and also ringdown transients occurring occasionally, typically due to servo-controls 
instabilities or mechanical relaxation in suspension system etc. While no precise 
model can be given for this noise until the detector is completed and fully tested, 
matched filtering techniques cannot be used to locate/remove these noisy signals. 

We propose a denoising method based on adaptive linear prediction techniques 
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which does not require any precise a priori information about the noise character- 
istics. Although our method does not pretend to optimality, we believe that its 
simplicity makes it useful for data preparation and for the understanding of the 
first data. 

In the following, we present the structure of the proposed algorithm and some re- 
sults obtained with the data from the Caltech 40 meter proto-type interferometer a. 
For a more detailed presentation, we refer the reader to □. 

2. Adaptive linear prediction 

The idea is to predict the current signal sample x k with a collection of past samples 
Xk = (xk—d-n, n = 0, 1, . . . , N — 1)*, the delay d > 1 being fixed arbitrarily. This is 
possible, only if the target sample shares enough information with (i.e., is sufficiently 
correlated to) the previous ones. In other words, the only predictable part of the sig- 
nal is the one whose correlation length is sufficiently large (i.e., long-term sinusoids 
or ringdowns). On the other hand, the broad band noise cannot be predicted, as it is 
not possible to guess the next value in this way. The prediction y k of x k is obtained 
through a linear combining of these data samples weighed by the corresponding 
coefficients w n , forming the tap-weight vector W — (w n ,n = 0, 1, ...,N— 1)', 
therefore leading to y k — W t X k . 

The optimal tap-weight vector W* which leads to the smallest prediction error 
efc = yk — Xk in the mean square sense can be proved to minimize a convex cost 
function. This minimization can be done using an approach similar to the steepest 
descent method with the following evolution equation for the tap-weights referred 
to as Least Mean Square or LMS algorithm i : 

W k+1 =W k + f ie k X kl (1) 

where the step gain parameter fi is an adjustable parameter. The tap-weight coef- 
ficients are renewed iteratively so that to converge and stabilize in a neighborhood 
of W* whose size is defined by /z. 

Once the filter has converged, we reject the predicted part of the signal (i.e., y k ) 
corresponding to the long-term sinusoids and the ringdown noise and we send the 
rest of the signal (i.e., e k ) for further analysis for detection. 

This LMS based prediction method is referred to as adaptive line enhancer 
(ALE) In this context, the term "adaptive" has two different meanings. First, 
it means that it will auto-adjust to reach for the best setup for a problem which 
is not initially precisely defined. Second, it is also able to follow changes in the 
characteristics of the data being processed in case they occur. 

Convergence time, frequency tracking ability and frequency resolution are con- 
trolled by the three adjustable parameters : the number of tap-weight coefficients 
N, the step gain parameter [i and the prediction depth d. One can take advantage 
of certain settings of the ALE to select a family of signals intead of another. 
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3. The ALE in practice 

Structure of the algorithm — We decompose the frequency axis in p disjoint 
frequency subbands of the same size. In each subbands, we apply twice the ALE 
with different sets of parameters. In the first step, the adjustable parameters are 
tuned to best remove long-term sinusoidal components of the noise ;_ whereas in 
the second one, the target is the short-time oscillatory transients (see □ for details 
about parameter adjustment). 

The ALE needs to be applied only in the parts of the signal which appears 
non-Gaussian. Some refinements are adjunct for this purpose : the acceptation 
(dismissal) of the first cleaning step relies on the detection of a long-term sinusoids 
of sufficient amplitude in the data. In the second step, the ALE is applied only 
if the filtered output ul deviates from Gaussianity. Details about these additional 
vetos can be found in □. 

Results on Caltech 40m proto-type data — We have applied the algorithm 
to the Caltech 40meter proto-type data taken in October 1994. Figure |l| illustrates 
how the algorithm is operating in the fifth frequency subband (from 617 Hz to 771 
Hz) among the p = 32 ones being processed. Figure || shows comparisons between 
the power spectra and histograms of the signal before and after denoising. 

4. Concluding remarks 

The originality of the proposed approach lies in the fact that it is possible to treat 
oscillatory transients. However, it remains that a comparison of the performances 
achieved by this algorithm on noise features of longer duration with other existing 
methods need to be done. Finally, although undertaken in □, this method suffers 
from the lack of a complete statistical characterization. 
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ing us the Caltech 40meter proto-type data. 



References 



1. J. D. E. Creighton. Listening for ringing black holes, jr-qc 9712044, 1997. 

2. B. Allen et al. Observational limit on gravitational waves from binary neutron star in 
the Galaxy. Phys. Rev. Lett, 83(8):1498-1501, 1999. 

3. A. Abramovici et al. Improved sensitivity in a gravitational wave interferometer and 
implication for LIGO. Phys. Letters A, 218:157-163, 1996. 

4. E. Chassande-Mottin and S. Dhurandhar. Adaptive filtering techniques for the removal 
of noisy oscillatory transients from interferometric data. In preparation, 2000. 

5. B. Widrow and S. D. Stearns. Adaptive Signal Processing. Prentice Hall, Englewoods 
Cliffs, 1984. 

6. The Caltech signal has been downloaded and calibrated with the GRASP package 
(version 1.9.3, http://www.lsc-group.phys.uwm.edu). 



4 Adaptive filtering techniques for interferometric data . 



(a) input signal, band #5 
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Fig. 1. Illustration of the denoising procedure «i Caltech proto-type data. In the 

subband #5 (between 617 Hz and 771 Hz), the signal Q (the data were taken on the October, 
14th 1994, frame #2) contains two power line harmonics (at 660 Hz and 720 Hz), as we see on 
the spectrogram (a) (this is a time-frequency representation of the signal energy. Dark regions 
are associated to large values of the energy density). We apply the ALE a first time to suppress 
long-term components (see spectrogram (b)) and a second run (c) eliminates artefacts of shorter 
duration (such as fast fluctuations in the harmonic envelope). 
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(b) histogram, input signal 
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Fig. 2. Comparison between power spectra and histograms of input/output signals. 

The figure depicts power spectra (left column) and histograms (right column) of the Caltech 40 
meter signal (top row) in the operating frequency band, between 300 Hz and 1kHz and the same 
signal after denoising (bottom row). The histograms are displayed in a graph with special axes 
where a Gaussian bell curve should appeared as a straight line. 



